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Abstract 


In this work, a design of transimpedance amplifier (TIA) for cryogenic scanning 
tunneling microscope (CryoSTM) is proposed. TIA with the tip-sample component 
in CryoSTM is called as CryoSTM-TIA. With transimpedance gain of 1 GQ, the 
bandwidth of the CryoSTM-TIA is larger than 300 kHz. The distinctive feature of 
the proposed CryoSTM-TIA is that its pre-amplifier is made of a single cryogenic high 
electron mobility transistor (HEMT), so the apparatus equivalent input noise current 
power spectral density at 100 kHz is lower than 4 (fA)?/Hz. In addition, “bias-cooling 
method” can be used to in-situ control the density of the frozen DX?” centers in the 
HEMT doping area, changing its structure to reduce the device noises. With this 
apparatus, fast scanning tunneling spectra measurements with high-energy-resolution 
are capable to be performed. And, it is capable to measure scanning tunneling shot 
noise spectra (STSNS) at the atomic scale for various quantum systems, even if the 
shot noise is very low. It provides a powerful tool to investigate novel quantum states 
by measuring STSNS, such as detecting the existence of Majorana bound states in the 
topological quantum systems. 


1 Introduction 


High performance transimpedance amplifier (TIA) for cryogenic scanning tunneling mi- 
croscope (CryoSTM) is a key element [ij 2]. Recently, in Ref.[3], it was proposed a 
design of TIA for CryoSTM. TIA with the tip-sample component in CryoSTM is called 
as CryoSTM-TIA. In that design, the CryoSTM-TIA has a transimpedance gain of 1 
GQ, a bandwidth of more than 300 kHz, and its equivalent input noise current power 
spectral density (PSD) that characterizes the circuit inherent noise is 21 (fA)?/Hz at 100 
kHz. The CryoSTM-TIA is capable of measuring the scanning tunneling shot noise spec- 
tra (STSNS) of various quantum systems at atomic resolution, even if the measured noise 
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current PSD is as small as a few (fA)?/Hz. In Ref.[3], it was illustrated how the CryoSTM- 
TIA is used to measure STSNS of quantum systems to investigate novel quantum states, 
such as detecting Majorana bound states (MBSs) in iron-based superconductors [4]. In 
the CryoSTM-TIA, cascade the pre-amplifier (Pre-Amp) and the post-amplifier (Post- 
Amp) to form an operational amplifier (OPA), and denoted as Macro-OPA. Pre-Amp in 
Macro-Amp is a differential amplifier made of a pair of high electron mobility transistors 
(HEMTs), where the HEMT is a kind of cryogenic GaAs MESFET [5} [6]. However, it is 
quite difficult to select cryogenic GaAs MESFETs with identical performances for pairing 
[2]. 

In this work, a design for the CryoSTM-TIA with a single HEMT in Pre-Amp is 
proposed. In the CryoSTM-TIA, cascade Pre-Amp and Post-Amp to form an inverting- 
amplifier (Inv-Amp), which is not OPA. The CryoSTM-TIA still has transimpedance gain 
of 1GQ and bandwidth more than 300 kHz. However, its equivalent input noise current 
PSD is only 3 (fA)?/Hz at 100 kHz, i.e. 1/7 of that in Ref.[3], since a noisy HEMT is 
reduced. Therefore, the lower tunneling shot noise of quantum systems can be measured 
at the atomic scale with this apparatus. Furthermore, only a single HEMT is used in 
Pre-Amp, avoiding the difficulty of matching two identical HEMTs. And, “bias-cooling 
method” can be used to in-situ reduce the inherent noise of the HEMT. With these 
advantages, this apparatus will be a powerful tool to investigate novel quantum states in 
various quantum systems, in the wider range of applications than that in Ref.|3], such as 
high-T, superconductors [7| [8], topological superconductors [4], MBSs [4 9] [12], 
Andreev reflection [14], and Kondo effect [15], etc. 


2 Circuit of the proposed CryoSTM-TIA 


The circuit of the proposed CryoSTM-TIA is shown in Figfi| It consists of several com- 
ponents: the single HEMT amplifier part of Pre-Amp shown in dashed box (al) of Figf1| 
the constant-current source part of Pre-Amp in dashed box (a2), the post-amplifier (Post- 
Amp) shown in dashed box (b), the compensated feedback network shown in dashed box 
(c), and the signal source circuit shown in dashed box (d). The two stage amplifier made of 
Pre-Amp and Post-Amp is called as inverting-amplifier (Inv-Amp). Inv-Amp is connected 
with the feedback network to form TIA. The input G of TIA is connected to the signal 
source circuit, to form the CryoSTM-TIA. The components placed in the cryogenic zone 
are shown in the dotted box. The parameters of all components of the CryoSTM-TIA 
circuit are listed in Table 


2.1 Design of Pre-Amp 


As shown in Fig[i] only a single cryogenic HEMT is used in Pre-Amp. The HEMT in this 
work is CNRS-HEMT (denoted as C-H) [5] [6] developed by CNRS/LPN in France with 
the excellent cryogenic and noise characteristics. It is capable of operating under 0.5 K 
with only 0.1 mW on the ideal operating point of “Va, = 100 mV and Ja, = 1 mA”. Its 
parameters are listed in Table j! where ex is its equivalent input noise voltage PSD and 
ie, is its equivalent input noise current PSD [6]. 

The single HEMT amplifier part of Pre-Amp is shown in dashed box (a1) in Fig/1] The 
source of C-H is grounded by a resistor R, and a small variable resistor Rg. Rg = Rs + Rs1, 


Figure 1: Circuit of the proposed CryoSTM-TIA. Single HEMT amplifier part of Pre-Amp 
is shown in dashed box (al), constant-current source part of Pre-Amp in dashed box (a2), 
Post-Amp in dashed box (b), compensated feedback network in dashed box (c), and signal 
source circuit in dashed box (d). The components placed in the cryogenic zone are shown 
in the dotted box. The parameters of all components of CryoSTM-TIA circuit are listed 
in Table 


and a capacitor Cs of 0.1 mF is in parallel with Rg. Drain O1 of C-H is connected to the 
load resistor R,, and Lj as the other end of R; is connected to output P of the constant- 
current source through a small resistor Rpi. The resistor Rp is grounded through the 
variable resistor Rnı. Ry = Ry + Ry. The other end of Rp, denoted as Og, is connected 
to the resistors Rg, and Lg as the other end of Rə is connected to P through a small 
variable resistor Rp2. Ry = R2 = Ry. Lı and Lz are grounded by the capacitors C1 and 
C2 respectively, and Cy = C2 = Cy, =10 pF. Ca is the input capacitance of Pre-Amp. 
C-H, Rs, Rn, R1, R2, and Rp) are placed in the cryogenic zone, shown in the dotted box. 
C'i is the capacitance of the cable that connects gate G of C-H and the tip in the CryoSTM. 
C-H is as close to the tip as possible to reduce Ct, which can be less than 0.5 pF [18]. The 
operating point adjustment for Pre-Amp is shown in Sect [4] The Pre-Amp parameters 
are shown in Table [I] 

For Pre-Amp, when AC voltage V is applied to input G, the AC voltage difference 
between O1 and Og is denoted as Vop- Ayp = Vop / V is the voltage gain of Pre-Amp. Since 
the gain-bandwidth-product of the CNRS-HEMT is gm /[27(Cgs+Cga)] ~ 1 GHz, the band- 
width of Pre-Amp is more than 30 MHz. In max{gm/(27Cg), 1/(27RgCg), 1/(2rRsC2)} << 


Table 1: Parameters of all components of CryoSTM-TIA 


circuit 
CNRS-HEMTs C-H 
Gate-source resistance Ra >10 TQ 
Transconductance gm 40 mS 
Channel conductance ga 1 mS 
Gate-source capacitance C's 5 pF 
Gate—drain capacitance Cga 1 pF 
Drain-source voltage Vas 100 mV 
Drain-source current Ias 1 mA 
e2 ((nV)?/Hz) 10 kHz 0.25 
100 kHz 0.07 
i} ((fA)? /Hz) 10 kHz 0.1 
100 kHz 1 
Pre-Amp 
Rs 100+5 Cs 0.1 mF 
Ruy 200 +5 2 RL 1 kQ 
Rpı 10 Q Rp2 10+10 Q 
C1& C2 10 uF 
Tı & T2 BJT BFT93 [16] R» 347+10 
Ry 20 kQ Vpp +12 V 
Post-Amp with THS4021 as Rear-OPA 
aao 97.5 dB fo 14.5 kHz 
Ca 1.5 pF Ra 1 MQ 
Re 330 kQ Ĉi 100 pF 
Supply voltages Vrp, Van +15, -15 V 
Feedback network 
Rp 1 GQ CF ~ 0.3 pF 
Rk 100 kQ Ck ~ 0.2 pF 
Ce 3 nF 
Signal source circuit 
Ry >1 MQ Ci ~ 0.5 pF 
Note: + indicates the variable resistance range. Without specifi- 


cation, the default value after + is 0. 


f <3 MHz (i.e. 112 Hz << f < 3 MHz with the parameters listed in Table (1), 
Ayp x —9m Ra, (2.1) 


where Ra = RL/(1 + Ruga). i.e., Avp ~ —20 in 112 Hz<< f < 3 MHz. The input 
capacitance of Pre-Amp is 
Ca = Cgs + (1 = Ayp)Cga- (2.2) 


i.e., Ca ~ 26 pF in 112 Hz<< f < 3 MHz. The input resistance of Pre-Amp Ra is the 
gate-source resistance of C-H, and Ra > 10 TQ, so it can be considered as infinity. 

The constant-current source part of Pre-Amp shown in the dashed box (a2) in Fig/1] 
is the same as in Ref.[3]. For a given resistance Rp, there is almost no fluctuation for the 
current Isour generated by the constant-current source, even though the voltage Vp, of 
the positive voltage source fluctuates greatly, which ensures the stability of the operating 
point for C-H. 


2.2 Design of Post-Amp and composition of Inv-Amp 


In Figf1] the Post-Amp circuit in the dashed box (b) is the same in structure as that in 
Ref. [B]. There is a commercial operational amplifier (OPA) in the circuit, called as Rear- 
OPA, which is the OPA with high gain-bandwidth-product, such as THS4021, OPA657, 
and LMH6624, etc. Ra and Ca are the input resistance and capacitance of Rear-OPA 
respectively. The feedback resistor Rs is connected to output O of Rear-OPA and its 
inverting input. Rr; is 330 kQ, and it placed in the cryogenic zone. Post-Amp also contains 
two cables that connect the non-inverting input and inverting input of the Rear-OPA to 
outputs O; and Og of Pre-Amp respectively. The capacitance of the two cables is Ci; and 
Cig respectively, and they may vary from 50 pF to 150 pF. In this work, Ci, = Cig = Ci = 
100 pF is assumed. The Rear-OPA in this work is THS4021 [17], and Ci >> Ca. The 
Post-Amp parameters are listed in Table [i] 

Cascade Pre-Amp and Post-Amp to form Inv-Amp. For the AC signal, the voltage 
gain of Inv-Amp is aa = aa(f). aa(f) can be expressed as 


QA = Ayp Ayr. (2.3) 


By the nodal analysis method, aa can be obtained. And then, with Ayp expressed by 


Eq. (2.1), Avr can be obtained by Eq. (2.3). In max{gm/(27Cs), 1/(27RgCs), 1/(27RgC2)} << 
f <3 MHz (ie. 112 Hz << f < 3 MHz with the parameters listed in Table i), 


Re 1 + j2r f RyuLCi 1 
Ru 1+j2rfRaCi 1+ 4 jon fe’ 


aa RyL 


Ayr ~ (2.4) 


where Ry, = RyRu/(Ru + RL) and aa is the voltage gain of the Rear-OPA. In this 
work, Rp >> Ry, and Ra >> Ry. And, a, can be approximately expressed as aa = 
aao/(1 + 9f/fp) in (0, 40 MHz] for THS4021. 
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Figure 2: The Inv-Amp voltage gain aa(f). The solid curves are the curves of |aa(f)|ap, 
and the dashed curves are the curves of Z(aa(f)). The red curves are the TINA-TI 
simulation results, the green curves are the calculated results with the nodal analysis 
method, and the black curves are the calculated results by Eqs. (2.1), (2.3), and in 
Article. 


With the parameters in Table[1} the Inv-Amp performances can be simulated by TINA- 
TI [19]. Fig[2| show the curves of |aa(f)|ap and Z (aa(f)) simulated by TINA-TI. By the 


nodal analysis method, aa(f) is calculated by the equations established with all compo- 
nents listed in Table |1| [20], which is also shown in Fig || The curves of jaa(f)|ap and 
Z (aa(f)) calculated by the nodal analysis method are identical with the simulated ones 
respectively. |aa(f)|ap and Z(aa(f)) are also calculated by Eqs. (2.1), (2.3), and (2.4), 
and the calculated curves in Fig [2] are almost identical with the simulated ones respectively 
in [1 kHz, 10 MHz], which verifies the correctness of Eqs. (2.1), (2.3), and (2.4). 

A 1 meter cable with characteristic impedance of 50 Q has a distributed capacitance of 
100 pF and a distributed inductance of 250 nH, and its distributed LC structure results in 
several resonant gain peaks in [80 MHz, 1 GHz] [21]. In Inv-Amp, Post-Amp and Pre-Amp 
are connected by the cables, and the resonant gain peaks may cause the self-oscillations, as 
some unexpected electromagnetic couplings may exist among the outputs and the inputs 
of the amplifiers. The simulations show that, by adding two low-pass filters (with the 
upper cut-off frequency of 3 to 10 MHz) at the Rear-OPA inputs or at the ends of the 
cables (depending on the type of low-pass filters being used), these resonant gain peaks 
can be greatly reduced, thus avoiding the self-excited oscillation. 


2.3 Frequency compensation of feedback loop 


In order to increase the bandwidth of the CryoSTM-TIA, for the high feedback resistor 
Rp with parasitic capacitance Cp, frequency compensation must be used in the feedback 
loop [I 22]. In Fig{i} the compensated feedback network is shown in the dashed box (c). 
Rp is placed as near as possible to C-H for reducing the capacitance between the gate 
G of C-H and Ground. Taking C, equalto kCp, where k is above 103, adjust Rẹ equalto 
Rp/k, realizing RgCe = RpCp. The output voltage of TIA as V, generates the current Ip 
flowing to input G of TIA, so 


Ze f) = 2 x oi ae N Bp 
2 İp 14+j2nfR.C, 1+ 52nfR,C,’ 


where Ci is the parasitic capacitance of Ry 22]. Zp(f) can be considered as the 
impedance of the feedback network. In Ref.{22], it has been achieved in experiments to 
broaden the bandwidth of the feedback network with the very high feedback resistor Rp 
of 10 GQ to MHz. In (0, 1 MHz], with the parameters listed in Table |1| |Zp(f)| ~ 
Rp/|1 + j27fRC| > Rp /1.008 and |Zp(f)| < Rp, so it can be considered that Zp(f) is 
equal to Rp. 

In order to verify the above frequency compensation method, the following experiment 
is done. A resistor with 500 MQ(+25 ppm/°C) at room temperature is chosen as Rp, which 
is 1.16 GO at 4.2 K. Rx is a 380 kQ resistor in series with a0 ~ 20 kQ potentiometer, and Ce 
is a COG ceramic capacitor of 10 nF. In order to easily adjust the frequency compensation, 
a 3 pF capacitor is paralleled with the Rp, so that the Cf is about 3.2 ~ 3.4 pF (estimated 
as 3.34 pF). Fig.3 shows the experimental results after frequency compensation, in which 
|Zp| = Rp/V2 at 905 kHz, and zp as the phase of Zp is —31.5° at 905 kHz. 


2.4 Circuit stability of the proposed CryoSTM-TIA 


In Figfi] the signal source circuit is shown in the dashed box (d), and its parameters are 
shown in Table The differential resistance of the tip-sample tunnel junction (TJ) is 
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Figure 3: At 4.2 K, Rp is 1.16 GQ and Cf is estimated as 3.36 pF. As Ry = 390 kQ and 
Ce = 10 nF, the black curve shows that the bandwidth of |Zp|/Rp is broadened to 905 
kHz, and the red curve shows zp as the phase of Zp is —31.5° at 905 kHz. 


Ry, which is limited to no less than 107? Rp in this design. The capacitance of TJ is Cy, 
which is in parallel with Ry. And, Cy is estimated as several fF [3]. C = Ca + C,+Cy in 
this work. Cy is at least two orders of magnitude less than Ca + C1, so it can be ignored 
in C and C ~ Ca + Cy. The DC bias & modulated signal voltage source is denoted as 
BMS, which provides the DC bias V; and sinusoidal modulated signal voltage V; for the 
CryoSTM-TIA. In the following simulation, C{ is always taken as 0.5 pF. TIA connects the 
signal source circuit to form the CryoSTM-TIA. According to the parameters in Table [I] 
the performances of the CryoSTM-TIA can be simulated by TINA-TI. 
The loop gain Ty, of the proposed CryoSTM-TIA [23] is 


TL(f) = aa (PEC) = -aa (FVU; 
in which 8(f) is the feedback factor, and its reciprocal is, 


1/B(f) ~ 1 + Z[1/R + 1/Ra + j2rf (Ca + C1)]. (2.5) 


By Eq. (2-5), |1/8(f)| and Z (1/6(f)) can be calculated. Fig[4(a) shows the calculated 
results of |1/6(f)|ag and < (1/G(f)) with Ry = +00, and Fig/4{b) shows those with 
Ry = 1 MQ. |aalap and Z (aa) are also shown in Fig] Both figures show |TL]aB = 
nalas =A as < 10 dB im PS 512 Me end 20.) = Zo) — 24/2) — 
Z (aa) — 180° — Z (1/8) > —138° in f < 512 kHz. Therefore, the CryoSTM-TIA is stable 
with gain margin more than 10 dB and phase margin more than 42°. 


2.5 Voltage gain and transimpedance gain of the proposed CryoSTM- 
TIA 


With the compensated feedback network mentioned in Sect [2.3] it can be considered that 
Zp is equal to Rp in (0, 1 MHz]. Considering the TJ capacitance Cy, the TJ impedance 
should be Zj = Ry/(1+ j2afRyCy). And, C ~ Ca + Cj. As the AC input voltage V; is 
applied by BMS, the output voltage of the CryoSTM-TIA is Vo, and the voltage gain of 
the CryoSTM-TIA is Ay = Vo/V;. In (0, 1 MHz], by the nodal analysis method, Ay is 
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Figure 4: TINA-TI simulation results for the voltage gain Inv-Amp a4, and 1/6(f) (a) with 
Ry = +00, (b) with Ry = 1 MQ. Both figures show |Tylap = |aalap — |1/6(f) lap < —10 
dB in f > 512 kHz and Z (TŁ) = Z(—a@a) — < (1/8) = Z (aa) — 180° — Z (1/8) > —138° 
in f < 512 kHz. Hence, the CryoSTM-TIA is stable enough. 


Setting V; = 0 and applying a sinusoidal current source I, in parallel with TJ, the 
output voltage V, is generated at the output of the CryoSTM-TIA. A; = Vo/J; is called 
as the transimpedance gain of the CryoSTM-TIA. In (0, 1 MHz], Aj is 


Rr 
Aj ood 1 1 Rp Rp 2r f Ryp (Ca +C) i (2.7) 
aA aa Ry aa RA J an 


Disconnecting TIA with the signal source circuit, and applying a sinusoidal current 
source Jip into the input of TIA, the output voltage Vor is generated at the output of 
TIA. Ait = Vor/Jir is called as the transimpedance gain of TIA. In (0, 1 MHz], Air is 


Rp 
Ait = (2.8) 
1 R . RrC,? 
1- as a RA — j2n f T 


In Fig)5| |Air/Relap and Z (Air/Rp) calculated by Eq. (2.8) with the parameters in 
Table[l]is identical with their TINA-TI simulation results in [1 kHz, 1 MHz], which verifies 
the correctness of Eq. (2.8) for Air. The upper cut-off frequency of Ajr(f)/Rp, i.e. its —3 
dB frequency is far = 320 kHz. Comparing Eq. (2.7) with Eq. (2.8), as Ca + Ci + Cy & 
Ca + Ci ~ Ca and Ry > 107*Rp, the upper cut-off frequency of the CryoSTM-TIA is 
approximately equal to fpr. 


AB m r rre errr e a ere aera a NEG 


— Mir'Rrlas 
—— Calculated |4;/Rplag 
=-= L(A /Rp) 

~~~ Calculated 2(4;,/R,) 


' 
n 
4-7 


-10F 


B rsrnnul sorral TT ssl oarl eT TTY 


10° 10° 10° 10' 10° 10° 10° 10° 10 


J (Hz) 


Figure 5: TINA-TI simulation results for Air/Rp, Air is the transimpedance gain of TIA 
in the proposed CryoSTM-TIA. |Air(far)/Relag = —3 dB and Z (Ait(far)/Rr) = 99° 
at for = 320 kHz. 


2.6 ‘Transient response of the proposed CryoSTM-TIA 


The simulation results for the transient response of the proposed CryoSTM-TIA are shown 
in Supplemental file 2 [24]. For the CryoSTM-TIA, the time taken from adding the input 
step signal voltage to the output response stably within a certain error is called transient 
response time tr. For the simulations, a resistor with a constant resistance Ro is instead 
of the tip-sample junction. When Ro > 30 MQ, for 0.1% error, tp < 5 us. When 5 
MQ < Ro < 30 MQ, for 1% error, ty < 5 us. When 1 MQ < Ro < 5 MQ, for 3% error, 
tr < 5 ps. 


3 Inherent noise of the proposed CryoSTM-TIA 


For the circuit of the proposed CryoSTM-TIA shown in Fig[i] the differential equivalent 
circuit with all noise sources is used to calculate its equivalent input noise. The details for 
the noise calculations are shown in Supplemental file 3 [28]. 


3.1 Equivalent input voltage noise and equivalent input current noise of 
Inv-Amp 


The equivalent input noise voltage and equivalent input noise current of C-H are denoted 
as ey and iy respectively. The resistors R1, Ro, and Rg are in the cryogenic zone of 4.2 K, 
and their noises in f > 1 kHz are thermal noise, which can be neglected [8]. The thermal 
noise voltage of the resistor Ry is ery. The equivalent input noise voltage and equivalent 
input noise current of the Rear-OPA are denoted as ea and ia respectively. These noise 
sources are independent. The equivalent input noise voltage and equivalent input noise 
current of Inv-Amp are denoted as ea and ia respectively. By the nodal analysis method 
and Wiener-Sinchin theorem, ignoring the minor terms, it is obtained that 


> =, hae e ( i2 
2 2 HL “RH a a 
e =e24 + +(1+—) 3, (3.1) 
A H Ae A’, A?p Ra g2, 


> a Cl (Ritz a R "E 
A = a ae (Qn f)? A ( HL eby { a) } (Qn f)? (Ce + Cga + Fen) g2? (3.2) 
H m 


RoR — — 
eaii =Cine,)” = jonf & *( e a) 
Aa 


(3.3) 


NY 
AN 
— 

+ 
z| y 

ai: 
NY 
repre) 


— j2r f (Cp + Cga + Pat oy 
Ra 

According to the parameters listed in Table |1| for Inv-Amp in this work, BA = 0.27 
(nV)?/Hz and a = 0.13 (fA)?/Hz at f = 10 kHz, and e = 0.08 (nV)? /Hz and A = 2T 
(£A)? /Hz at f = 100 kHz. For the Macro-OPA in Ref. [3], e2 = 0.5 (nV)? /Hz and A = 0.8 
(£A)? /Hz at f = 10 kHz, and BA = 0.14 (nV)?/Hz and B = 20 (fA)? /Hz at f = 100 kHz. 
The resistor Ry in Inv-Amp in this work replaces the transistor H2 in the Macro-OPA in 
Ref. [B]. The noise generated by Ry is 2 orders of magnitude smaller than that generated 
by H2 in the Macro-OPA in Ref.[8], i.e. Coe << Are. Therefore, the equivalent input 


noises of Inv-Amp in this work are much lower than those of the Macro-OPA in Ref.[3]. 
In f > 10 kHz, e2 = 2.25 (nV)?/Hz and i? = 4 (pA)?/Hz [IZ]. In Eqs.(8-1), 
(3.2), and (3.3), eby RAL / R3 is one order of magnitude smaller than e. In Eq. (8-1), 
(1+ Rut/Ra)? iZ/g2 is 2 orders of magnitude smaller than e2. In Eq. (3.2), (Ces + Coa + Calan Ray 12/92 
is one order of magnitude smaller than CR e2/A2 p: Ihn Eq. (3.3), (Cys + Cga + Ca Rut/Ra) (1+ Ru/Ra) 12/92, 


is one order of magnitude smaller than Cae? /A?p. Further ignoring the minor terms in 


e3 = ef + 2/Atp, (3.4) 
= i+ ~ CREA, (3.5) 
eata = (iach) = —j2n fC ae? /App- (3.6) 


3.2 Equivalent input current noise of the proposed CryoSTM-TIA 
The equivalent input noise current PSD of the proposed CryoSTM-TIA is obtained as 


i2 =i + 4kgT/Rr + [1/R} +1/R2 + (2r f) CR] e% 
+ (1/Ry + j2nfCy) (eath + €/Rr) 
+ (1/Ry — j2n fC) (Treh + sh | (3.7) 
where Ciy = C1 + Cy [3] [25] [28]. Putting Eqs. (3.4| (3.4), 6-5 , and (3.6) into Eq. (3.7), 
iZ =i, + 4kgT/Rp + (1/Ry +1/Rp)? e g =- 


+ (Qn f)? (c 2 e2 + C2e2/A2 Fl (3.8) 


2 =i + 4kpT/Rp + (1/Ry + 1/Rp)} e + iby, (3.9) 
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where i2,, = (2r f)?’ (Chei + 2CaCyye2/A2p). 

For the proposed CryoSTM-TIA, Rp = 1 GQ, Ca = 26 pF, Cy = 0.5 pF, and Cj = 10 
fF. Rp and TJ are in the cryogenic zone at 4.2 K. As Ry = 1 MQ, i?, and its four 
components are listed in Table [2] The noise components of the CryoSTM-TIA in Ref.[3] 
are also shown in Table The CryoSTM-TIA proposed in this work has an equivalent 
input noise current PSD of 0.6 (fA)?/Hz at 10 kHz and 3 (fA)?/Hz at 100 kHz, which is 


much lower than that in Ref.|3}. 


Table 2: Noise components of CryoSTM-TIAs 


Pre-Amp type Single HEMT Amp Pre-Amp in Ref. [3] 
f (kHz) 10 100 10 100 
e? ((nV)*/Hz) 0.27 0.08 0.5 0.14 

Unit for the following terms is (fA)?/Hz 

B 0.13 2.7 0.8 20 
4kpT/ Rr 0.2 0.2 0.2 0.2 
iZi 0.0013 0.06 0.03 0.7 
ERA 0.26 0.08 0.5 0.14 
i2, as total 0.6 3 1.5 21 


in 


For the Epi-wafer of the CNRS-HEMT, as Si is doped in Al,Ga,_,As, each dopant 
gives rise to two types of electronic states: a shallow and delocalized donor level associ- 
ated with the normal substitutional site configuration, called as shallow donor state; a 
more localized acceptor level, known as DX~ center, arising from a lattice distortion. In 
Al,Ga ,_,As, as x > 0.22, DX~ as a bound state is more stable than the shallow donor 
state. And, x = 0.37 for CNRS-HEMT in this work. There is a large repulsion barrier 
Ecap in the transition from the shallow donor to DX~, so DX? is frozen as T < 120 K 
. 

According to the work on the noise mechanism of the CNRS-HEMT [5], “bias-cooling 
method” can be used to control the density of the frozen DX~ centers in the doping area in 
the HEMT. This unique method, by which the in-situ modification of the HEMT structure 
is realized, can be used to reduce the noise generated by the HEMT. By cooling the HEMT 
with the positive gate-source voltage, more electrons in the HEMT are frozen into the DX~ 
centers. At low temperatures, the absolute of the gate-source voltage required to maintain 
the HEMT at the ideal operating point is reduced, so that the gate leakage current is 
reduced, and thus the low frequency noises generated by the gate leakage current is reduced 
[5]. By the means, the inherent noise of the CryoSTM-TIA can be reduced. Moreover, 
if the HEMT used for the CryoSTM-TIA is not CNRS-HEMT, or its noise performances 
are not as good as those listed in Table [I] it can be remedied by “bias-cooling method”. 
This approach is very convenient for cases with only a single HEMT in Pre-Amp. 


4 Operating state adjustment and DC tunneling current 
measurements 


For the CryoSTM-TIA in Fig[i] along with the circuit parameters listed in Table |1} its 
operating state can be adjusted as follows. (1) Disconnect Pre-Amp from Post-Amp and 
ground the input of Pre-Amp, i.e. gate G of C-H. Adjust Rp2 to Rp2 = Rpi. Adjust the 
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constant-current source to generate current as Isour = 2 mA. Adjust Rsı and Rp1 to achieve 
C-H at the ideal operating point (Vas = 100 mV, Jg, = 1 mA) and the potentials at O1 and 
Oz equal. That is to say, Rs; and Ry; are adjusted to realize Ry = Rs + Vas/Ias = Rg +100 
Q. (2) Cascade Pre-Amp and Post-Amp to form Inv-Amp, and input G is still grounded. 
By then, the DC output voltage Vom of Inv-Amp is usually not 0. Vom is induced by the 
input offset voltage of the Rear-OPA [I7], its input bias currents and input offset current, 
and the common-mode DC voltages on the the inputs of the Rear-OPA. Adjust Ry, and 
Rp2 to achieve the current through Rp2 back to 1 mA and Vom = 0. Therefore, C-H is 
also at the ideal operating point. (3) Connect output O and input G of Inv-Amp with 
the feedback network and disconnect input G from ground, to form TIA. Input G is still 
not connected with the signal source circuit, so the potential of input G and output O of 
Inv-Amp is still 0, since the input resistance of the transistor C-H Ra can be considered 
as infinity. In the feedback network, Rp +R, ~ Rp = 1 GQ. (4) Connect the signal source 
circuit to TIA, to form the CryoSTM-TIA. 

When f — 0, aa > aao, and aao is the DC voltage gain of Inv-Amp. |aaolan © 
(gmRe/gaRu)ap ~ 96.4 dB, which is consistent with the simulation result |aao|an = 95.6 
dB shown in Fig 2 As the DC bias V; is applied by BMS, the DC resistance of TJ is 
R, and the potential at input G is Va, and the output voltage of the CryoSTM-TIA is 
Vo. Obviously, aaoVa = Vo, and (Vi — Va)/R = (Va — Vo)/Rr. The DC bias on TJ is 
V =V,— Vg, and the DC tunneling current I = (V; — Va )/R. I; as an approximate value 
of I is 

I, = —V,/ Re. (4.1) 


And, the relative error is obtained as 
Er = |I; — I|/|I| = 1/(1 — aao). 


As |aaolap © 96 dB, Er < 20 ppm, consistent with the simulation results [29]. Since 
aaoVa = Vo, the DC bias on TJ is 


V == Wein: (4.2) 


By Eq.(4.1), Vo = —Rels ~ —RpI = —RpV/R. By Eq.(4.2), Vi ~ V — RpV/(Raao). As 
the minimum of Ry is not less than 1 MQ, R > 1 MQ, so |Rp/(Raao| < 1/60. Thus, 
V & Vi. With the measured V,, I ~ J, and V & V; can be obtained. Therefore, the 
scanning tunneling current spectra I = I(V) (V €[V,, Vy) can be obtained, where Vg 
and Vy are the lower limit and upper limit voltages for the measurements respectively. 

It should be noted that Pre-Amp in this work is not a differential structure, so there 
is no function for it to counteract the influences of the temperature drift as that for the 
differential structure. The drift of Vom can be considered as the amplification of the drift 
of the Inv-Amp input offset voltage Vos, where Vom = a@aoVos. The estimated drift of 
|Vos| is about 0.5 mV/°C [30]. The power of the constant-current source is about 24 mW. 
The well-designed temperature Control System based on the TEC devices can be used 
to control the temperature fluctuations of the constant-current source within 0.0002 °C 
and those of Rear-OPA within 0.1 °C, so the fluctuations of Vog within 150 nV is 
guaranteed, i.e. the TIA output fluctuations within 150 nV [3]. 

In order to prove the correctness of our proposed method for designing TIA, the N- 
Channel JFET (SST4393) is chosen to replace CNRS-HEMT, and TIA is designed and 
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fabricated with the same circuit design method and the same circuit topology. SST4393 
can work at 77 K. We measured the gain, noise, and transient response of JFET-based 
Inv-Amp and TIA at room temperature and 77 K [32]. The experimental results are 
basically consistent with the theoretical calculations and simulations, which shows that 
the method of designing TIA is correct. In addition, the designed JFET-based TIA can 
also be used for STM at 77 K, with larger bandwidth and lower inherent noise, comparing 
with traditional TIAs for CryoSTM (e.g. FEMTO DE-DLPCA-200 TIA). 


5 Applications for the proposed CryoSTM-TIA in spectra 
measurements 


For most applications, the modulus of the transimpedance gain of the CryoSTM-TIA 
|A;(f)| should be measured firstly. Aj;(f) is expressed as Eq.(2.7). Since |aa(f)lap > 90 
dB in [1 kHz, 300 kHz] as shown in Fig(2| when Ry > 107? Rp, |Aj(f)| in [1 kHz, 300 kHz] 
can be approximately expressed as 


|Ai(f)| = Re/|1 — j2r fRp(Ca + Cr)/aa(f)]- (5.1) 


How to measure it has been described in Ref. [3]. 


5.1 Measurements of scanning tunneling differential conductance spec- 
tra by the proposed CryoSTM-TIA 


The differential conductance of TJ Gy = 1/Ry is the function of the voltage V applied 
to TJ. As the frequency f of the modulated signal voltage V; provided by BMS is low 
enough, such as f < 1 kHz, by Eq. (2.6), 
Rr 1 
Ay & . ; 
RyV) 1- NE AAA 
Therefore, Rj(V) ~ [1/Ay +1/aa(f)] Rr. When the measured |Ay| < 1000, Ry(V) = 
Rp/|Ay|, since |aa(f)|ap > 85 dB in f < 300 kHz as shown in Fig[2| With the measured 
|Ay| = |Vo|/|Vi|, the differential conductance spectra Gy(V) = 1/R3(V) (V €[Vi, Vaļ) can 
be obtained. 
Increasing the frequency f of the modulated signal can speed up the scanning tunneling 
differential conductance spectra measurements. In [1 kHz, 300 kHz], as Ry > 1 MQ, by 


Eq. (2.6), 


(AIS po * 
. [Z3(f)| |1 j2rfRe(Ca + Cr)/aa(f)| 

Therefore, 1/|Z3(f)| ~ |Av(f)|/|Ai(f)|, where |Aj(f)| is shown by Eq.(5.1). And, |A,(f)| 
and |Aj;(f)| can be measured. 1/|Z3(f)| = 1/R3 + (2nfCy)* can be obtained. Selecting 
two different frequencies fı and fz in [1 kHz, 300 kHz], |Z3(f1)| and |Z3(f2)| are obtained. 
Ry and Cy can be solved out from the measured |Z 3(f1)| and |Z3(f2)|. 

Since the inherent noise of the CryoSTM-TIA is very small, the amplitude of the 
modulated signal voltage V; can be very small (< 10 uV), so that the energy resolution 
for STS is much improved. 
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5.2 Measurements of scanning tunneling shot noise spectra by the pro- 
posed CryoSTM-TIA 


The measurement method for the tunneling shot noise spectra by the proposed CryoSTM- 
TIA is basically the same as that introduced by Ref.[3], and it is only briefly described 
here. 

Before approaching the tip to the sample in the CryoSTM, Ry can be considered as 
infinity and Cy as 0. In this case, the output noise voltage PSD of the CryoSTM-TIA 
Ssu(f) can be measured, and the equivalent input noise current PSD i?(f) is 


P(f) = Seu A (5.2) 
where |A;(f)| is obtained by Eq. (5.1). 


To measure the scanning tunneling shot noise spectra (STSNS) of a quantum system by 
the CryoSTM-TIA, the distance between the tip and sample is adjusted, and the interval 
(VL, Vu] is selected, so that the shot noise measurements at V € Dy = {V| V% < V < 
Vu,Gy3(V) < 1 uS} are needed to study the physical properties of the quantum system 
[3]. In Dy, the output noise voltage PSD Ssum(f,V) can be measured. The tunneling 
noise current PSD is denoted as S}(f,V). The equivalent input noise current PSD of the 
CryoSTM-TIA is the function of Ry, so it is the function of V and denoted as i? (f, V). 
For Ssum(f, V), 


Sml V) = (V) BV| LAA? (5.3) 
With Eq. (5.2) and (5.3), it is obtained that 


sirv) = See al) _ Rv) -RA (5.4) 


From Eq. (3.9), 


po = ol 
+ (Qf)? (az + ES . 
vP 


> m 


6(f,V) as the difference between i2 (f, V) and i?(f) is smaller than 0.3 (fA)?/Hz in [10 
kHz, 100 kHz) and 0.1 (fA)?/Hz in[100 kHz, 300 kHz]. 5(f,V) can be neglected compared 
with S1(f,V), so long as the minimum of G3(V) (V € Dy) is not too small. Therefore, 
Si(f, V) can be obtained as 


And i? ( f,V) is composed by the four parts as shown in Eq.(3.9). According to Table |2] 


Ssum V — Ssu 
Sf, V) = ae w. 


by the measured Ssum(f,V), Ssu(f), and |Ai(f)|. And then, the SNTNS Sis(I) = 2FelļI| 
can be extracted from S{(f,V). Here, J = I(V) (V € Dy), and I is the DC tunneling 
current as the bias V is applied. 

Compared with the CryoSTM-TIA in Ref.[3], the CryoSTM-TIA proposed in this 
work has the same bandwidth and transimpedance gain, but its inherent noise is much 
lower. The inherent noise of the CryoTM-TIA in this work is only 3 (fA)?/Hz at 100 kHz 


(5.5) 
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and its 6(f,V) < 0.1 (fA)?/Hz at 100 kHz, while that of the apparatus in Ref.[3] is 21 
(fA)?/Hz at 100 kHz and ô( f, V) < 0.28 (fA)?/Hz at 100 kHz. Therefore, as investigating 
novel quantum states of various quantum systems, the measurements performed with the 
apparatus proposed in this work is more accurate. 

For example, as investigating the existence of MBS in a magnetic flux vortex of the 
iron superconductor in CryoSTM [4], the tunnel junction resistance must be large enough, 
since the occurrence of incoherent Andreev reflection can be ruled out by weak tunnel 
coupling conditions [84]. In Ref.{4], the tunneling current I is quite low (see Fig.3(a) in 
Ref.{4] and Fig.S1(a) and (b) in its supplemental file). As the tunnel junction bias V is 
0.3 mV, J is only tens of pA, and the corresponding shot noise may only be only a few 
(fA)?/Hz. For the shot noise measurements in the system in Ref.[4], it is obvious that the 
CryoSTM-TIA proposed in this work is much more effective and accurate than that in 
Ref. [3]. 


6 Conclusion 


In this work, a design of transimpedance amplifier (TIA) for cryogenic scanning tunneling 
microscope (CryoSTM) is presented. TIA connected with the tip-sample component in 
CryoSTM is called as CryoSTM-TIA. The CryoSTM-TIA in this work has transimpedance 
gain of 1 GQ and bandwidth more than 300 kHz, and its equivalent input noise current 
PSD is only 3 (fA)?/Hz at 100 kHz. In this CryoSTM-TIA, a single CNRS-HEMT Pre- 
Amp is instead of the differential Pre-Amp consisting of a pair of CNRS-HEMTs in Ref. [3]. 
The difficulty of matching the identical HEMTs is avoided. And, due to the reduction of 
a noisy HEMT, the apparatus inherent noise is only 1/7 of that in Ref.[3]. Furthermore, 
the inherent noise of the single HEMT in the circuit can be in-situ reduced by “bias- 
cooling method”. With this apparatus, the fast high-energy-resolution scanning tunneling 
spectra measurements can be performed and the very low tunneling shot noise of quantum 
systems can be measured at the atomic scale. This apparatus can be applied to investigate 
novel physical properties of various quantum systems, such as detecting the existence of 
Majorana bound states in the topological quantum systems. 
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Supplemental file 1: Voltage gain of Inv-Amp 


Connect to Constant Current-source 


Figure sl-1 The Inv-Amp circuit with all components for Eqs.(s1.1-s1.8), and their values are 


shown in Table 1 in Article. 


With the nodal analysis method, the Inv-Amp voltage gain aa(f) is calculated by the 
following equations. All components for Eqs.(s1.1-s1.8) are shown in Fig.sl-1, and their values 


are shown in Table 1 in Article. 


(V, -V,)/Rp = (Vs Vig) [Ro + Pa (Ry + Ry) + J@CY, (s1.1) 
joCy(V.-V,)=8,,(V.-V)+25(.-V,)+(V,-V,)/R+jaCh, , 612 
(V,-V,,)/R, =(Vg—Vo) [Ror + JOCK, l (s1.3) 
(Ta) [Ra = (Von Ve [Roo + JOC Vy (s1.4) 
(V,-Vo)/ Ry + (Vn Vo) / Ryo =9 (s1.5) 
8n(¥.-V,)+8,(V,-V,)=V,/R, + joCV, + jaC,,(V,-V,) (s1.6) 
V.=a,(V,-V,) , (s1.7) 


a, =V VÝ, n (s1.8) 
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Figure s1-2 The Inv-Amp voltage gain aa(f). The solid curves are the curves of |aa(f)|as, and the 
dashed curves are the curves of Z(aa(f)). The red curves are the TINA-TI simulation results, the 
green curves are the calculated results with the nodal analysis method (i.e. Eqs.(sl.1-1.8)), and 


the black curves are the calculated results by Eqs.(2.1), (2.3), and (2.4) in Article. 


Fig. s1-2 shows the modulus of the Inv-Amp voltage gain |aa(f)|ap. The red solid 
curve is the TINA-TI simulation results, and the green solid curve is the calculated results 
with the nodal analysis method. They are identical with each other. The black solid curve 


is the calculated results by Eqs.(2.1), (2.3), and (2.4) in Article. In fe[1 kHz, 3 MHz], 


three curves are identical with one another. 


Fig. sl-2 also shows the argument of the Inv-Amp voltage gain Z(aa(f)). The red 
dashed curve is the TINA-TI simulation results, and the green dashed curve is the 
calculated results with the nodal analysis method. They are identical with each other. The 
black dashed curve is the calculated results by Eq.(2.1), (2.3), and (2.4) in Article. In fe[1 


kHz, 3 MHz], three curves are consistent with one another. 


Supplemental file 2: Voltage gain peaks for cables 
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Figure s2-1 The circuit scheme for measureing the voltage gain peaks for cables. 


A 1 meter cable with characteristic impedance of 50 Q has a distributed capacitance of 100 pF 
and a distributed inductance of 250 nH, which is shown as Fig. s2-1. In Fig. s2-1, R=1.5 kQ, the 
gain is Vo/Vin. The simulation results show its distributed LC structure results in several resonant 
gain peaks in [80 MHz, 1 GHz] shown as Fig. s2-2. We amplify Vo with Amplifier SA-220F5 
(gain=200, fr=200 MHz), and then the gain peaks for the cable are measured with an oscilloscope. 


Fig. s2-3 show the measured gain peaks for cable. 
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Figure s2-2 The simulation results of the gain peaks for a 1 meter cable with a distributed 


capacitance of 100 pF and a distributed inductance of 250 nH. 
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Figure s2-2 The measured results of the gain peaks for a 1 meter cable with a distributed 


capacitance of 100 pF and a distributed inductance of 250 nH. 


Supplemental file 3: Transient response of the proposed CryoSTM-TIA 
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Figure s3-1 A resistor with a constant resistance Ro is instead of the tip-sample junction for the 
simulation. TINA-TI simulation results for the transient response of the CryoSTM-TIA with the 
different Ro for the time interval of 10 us. The dashed curve is the step input signal Vi/Vins: for the 
different Ro. The solid curves are the output response Vo/Vost. Here, Vost and Vins are their values 
respectively at 300 ms after applying the step input signal Vi, which are shown in Fig.s2-2. 


Transient response time ¢,<5 us. (b) is the zoom-in for the output response Vo/Vos in (a). 
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Figure s3-2 A resistor with a constant resistance Ro is instead of the tip-sample junction for the 
simulation. TINA-TI simulation results for the transient response of the CryoSTM-TIA with 
different Ro for the time interval of 300 ms. (a) The step input signal Vi/Vins. (b) The output 
response Vo/Vost. Vo/Vost= 1 in t= 100 ms. (c) is the zoom-in figure of (b). 


A resistor with a constant resistance Ro is instead of the tip-sample junction for the 
simulation. Fig.s3-1 shows TINA-TI simulation results for the transient response of the 
CryoSTM-TIA with the different Ro for the time interval of 10 us. The dashed curve is 
the step input signal V;/Vins, and the solid curves are the output response Vo/Vos for the 
different Ro. Here, Vost and Vinst are their values respectively at 300 ms after applying the 
step input signal Vi. The error is definite as |Vo—Vost |/|Vost|. For the CryoSTM-TIA, the time 
taken from adding the input step signal voltage to the output response stably within a 
certain error is called transient response time f. As shown in Fig.s3-1 and Fig.s3- 
2(b)&(c), with Ro= 1 GQ and 30 MQ, for 0.1% error, <5 us. With Ro=5 MQ, for 1% 
error, 4< 5 us. And, With Ro= 1 MQ, for 2.5% error, tr<5 us. 


Supplemental file 4: Noise of the proposed CryoSTM-TIA 


S4.1 Noises of STM-TIA 
For Inv-Amp Inv as shown in Fig.s1l-1, the equivalent input noise voltage of the Inv-Amp is 


denoted as ea and its equivalent input noise current is ia, and their harmonic components of 


E 
frequency f are Ea and J, respectively. ) can be obtained by the nodal analysis method. By 


A 


e ei, EE Er 
Wiener-Khintchine theorem, &s 8!) can be obtained ron iy rr [S4R1, S4R2]. 
A A 


eA 2 
LG. ty A A 


The two matrix elements on the main diagonal are the equivalent input noise voltage PSD of the 


Inv-Amp e% and its equivalent input noise current PSD i, . The two matrix elements on the 
sub-diagonal are its equivalent input noise voltage-current PSD e,i, and equivalent input noise 


current-voltage PSD i,e, . 


The Inv-Amp is connected to the feedback resistor Rr to form a TIA. The equivalent input noise 
voltage of the TIA is denoted as er and its equivalent input noise current is it, and their harmonic 
components of f are Er and Ir respectively. The temperature of the feedback resistance Rr is T, and 


its noise voltage is er and its harmonic component of f is Er. 


(b) 


Figure s4-1 (a) TIA circuit with the input short-circuit containing the equivalent input noise voltage of 
Inv-Amp Inv ea and its equivalent input noise current ia, and the output noise voltage of eortv; (b) Noiseless 
TIA circuit with the equivalent input noise voltage of the TIA er as the input signal, and the output noise 


voltage of eotve ; the equivalency of the above two circuits means @oty = €otve - 
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For the TIA, the circuit containing all noise sources with the input short-circuit is shown as 
Fig.s4-1(a), and the output noise is eoty. The noiseless circuit with the equivalent input noise voltage 
of the TIA er as the input signal is shown as Fig.s4-1(b), and the output noise is @otve. For 
calculating the equivalent input noise voltage of the TIA, the equations are established on the 
equivalency of the above two circuits, i.e. eoTv = @otve. Therefore, by nodal analysis method, 


Ertr= Ea. 


Figure s4-2 (a) TIA circuit with the input open-circuit containing the equivalent input noise current of 
Inv-Amp Inv ea and its equivalent input noise current ia, and the output noise current of eori, (b) Noiseless 
TIA circuit with the equivalent input noise current of the TIA ir as the input signal, and the output noise 


voltage of eotie ; the equivalency of the above two circuits means @oti = eotie . 


For the TIA, the circuit containing all noise sources with the input open-circuit is shown as 
Fig.s4-2(a), and the output noise is is eori. The noiseless circuit with the equivalent input noise 
current of the TIA ir as the input signal is shown as Fig.s4-2(b), and the output noise is eotie. For 
calculating the equivalent input noise current of the TIA, the equations are established on the 
equivalency of the above two circuits, i.e. @ori = eotie. Therefore, by the nodal analysis method, 


Tp=Iy+(Eat EpyRe . 


lerle] 


Hence, 


* * 2 F 
From e by Wiener-Khintchine theorem, ka Cri can be obtained. The two 
LE, hh Le; i 


matrix elements on the main diagonal are the equivalent input noise voltage PSD of the TIA e and 


its equivalent input noise current PSD 722. The two matrix elements on the sub-diagonal are its 
q p i g 
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equivalent input noise voltage-current PSD e,j, and equivalent input noise current-voltage PSD 


ase, (s4.1) 
=i +4k T/R +e R, (84.2) 
ei (ie) ei teh /R (54.3) 


Here, the noise voltage PSD of Rr is e = 4k,TR, . 


The TIA is connected with the signal source circuit to form a STM-TIA. The equivalent input 


noise current PSD of aSTM-TIA i can be obtained [S4R1-S4R3] by 

in =i +(1/R +(2nf) Ci) +(Y/R, + J2mfC, Jeri +(1/R,—j2mfC ine, (84.4) 
Here, Cy = Cr+ Cj. Putting Eq.(s4.1), ($4.2), and (s4.3) into Eq.(s4.4), A is 

in =i, + 4k? /R, +e, (I/R? +1/R +(2nf) Cz) 


Sa a (s4.5) 
+(I/R, + j2mfCy (ends +e, /Re) + (YR, J2mfCy) ine +e / Re] 


S4.2 Noises of the proposed CryoSTM-TIA 


Figure s4-3 Circuit of the proposed CryoSTM-TIA. Inv-Amp of the proposed CryoSTM-TIA is shown in 


dotted box. 
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For the proposed CryoSTM-TIA shown as Fig.s4-3 (i.e. Fig.1 in Paper), its Inv-Amp is shown 
as Fig.s4-4. For the Inv-Amp, its A P, ei, and e can be calculated by the nodal analysis 
method and Wiener-Khintchine theorem. Putting the results into Eq.(s4.5), the equivalent input 


noise current PSD of the proposed CryoSTM-TIA A is obtained. 


Figure s4-4 Inv-Amp circuit of the proposed CryoSTM-TIA. G is the input and O is the output. The 
single-transistor circuit in Pre-Amp is shown in dashed box (a1), the constant current source part in Pre-Amp 


is in dashed box (a2), and Post-amp is in dashed box (b). 


$4.2.1 The equivalent input noise voltage and the equivalent input noise current of Inv-Amp 
For the Inv-Amp, the equivalent input noise voltage and equivalent input noise current of 
transistor C-H are denoted as ep and iy respectively, and their harmonic components of f is are 
denoted as Ey and Jy respectively. gm is the transconductance of C-H and ga is channel conductance 
of C-H. The noise voltage of the resistors Ru is denoted as ern, and the harmonic components of f 
are denoted as Ern. The noise voltage of the resistors Ri, R2, and the feedback resistance Rr are 
denoted as e11, er2 and er respectively, and the harmonic components of f are denoted as E,1, E12, 
and Er respectively. Here, Ri= Ro= Rv. All noises of Ru, Ri, R2, and Rf above 1 kHz are thermal 
noise. Ra = Ri//(1/ga) and Rut = Ru//Ri are denoted. The equivalent input noise voltage and 
equivalent input noise current of Rear-OPA are denoted as ea and ia respectively, and the harmonic 
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components of f are denoted as Ea and J, respectively. These noise sources are independent of each 
other. For the Inv-Amp consistent of the Pre-Amp and Post-Amp, its equivalent input noise voltage 
and equivalent input noise current as ea and ia, the corresponding harmonic components of f are Ea 


and J, respectively. 


(a) + 
D i tm 


RH 


+ 
Cas 
o> Coy 

= me © eu 
l Cr XO f S) 
l + + 
TO 

Gi 
ey C 
O) 


Figure s4-5 (a) Inv-Amp equivalent differential circuit with the input short-circuit containing all noise 
sources, and the output noise current of eav; (b) Noiseless Inv-Amp circuit with the equivalent input noise 
voltage of the Inv-Amp ea as the input signal, and the output noise voltage of eoave; the equivalency of the 


upper two circuits means @pay = @oAve . In (a) and (b), the triangle OPA is the Rear-OPA. 


The equivalent differential circuit of the Inv-Amp containing all noise sources with the input 
short-circuit is shown as Fig.s4-5(a), and the output noise is eoay. The noiseless circuit with the 


equivalent input noise voltage of the Inv-Amp ea as the input signal is shown as Fig.s4-5(b), and the 
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output noise is @oAve. For calculating the equivalent input noise voltage of the Inv-Amp, the 
equations are established on the equivalency of the upper two circuits, i.e. @oav = @oAve. Therefore, by 
the nodal analysis method, 
Ea EqtErRuZu/(Avp; Ru)-Evi/(gmRv)tE1L2Zur/(AvejRi) 

+E Zu/R¢-Ea/Avpj —(A/gmt+Zur/Avej)la , (s4.6) 
where, Avpj = 2m/[1/Ritgay(2nfCi], Zur = Ruv/[1+ j(2n/RuCi]. 


Figure s4-6 (a) Inv-Amp equivalent differential circuit with the input open-circuit containing all noise 
sources, and the output noise voltage of eoai; (b) Noiseless Inv-Amp circuit with the equivalent input noise 
current of the Macro-OPA ia, as the input signal, and the output noise voltage of @oaic; the equivalency of the 


upper two circuits means @oai= @oaie. In (a) and (b), the triangle OPA is the Rear-OPA. 


The equivalent differential circuit of the Inv-Amp containing all noise sources with the input 
open-circuit is shown as Fig.s4-6(a), and the output noise is eoai. The noiseless circuit with the 
equivalent input noise current of the Inv-Amp i, as the input signal is shown as Fig.s4-6(b), and the 


output noise is eoaie. For calculating the equivalent input noise current of the Inv-Amp, the equations 
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are established on the equivalency of the upper two circuits, i.e. @oai= odie. Therefore, by the nodal 
analysis method, 
Ia = Tt j(2f)CajEruZuv/(Ayej Ru) -ATA Ces tC ga) ELER) + j(21f) CajEi2Zuv/(Avej Rv) 
+ j(2mf)\CajEZuv/(AvejRs) -JOSCE Avr j— OIC Cest Coa)/gim + CajZur/Avejla, (84.7) 
where, Caj = CystCga(1+Ayp)). 
Thus, 


a) (ela) Cree at benera iene, te 


1 Zak; 1 E, 18m +Z /4r 7 
(, i(2nf')C A en j(2af) (Ca +C, 1) RL N i 


(s4.8) 


For the measured frequency of fm < 300 kHz, 2rfCi <0.2 mS with Ci = 100 pF, while 1/Ra= 
1/Ri+ga=2 mS and 1/RuL = 1/Ri+1/Ru = 6 mS, so Avrj = gm/[1/Ri+gat7(20/)Ci] is approximately 
simplified to —Avp = gm/(1/Ri+ga) = gmRa, Zur to Rut/(gmRa), and Caj to Ca. Ignoring the small 


quantities, such as the thermal noise of R1, R2, and Re, 


ey =e +] (Ri. /Ra enn +e |/4. +(1+ Ry /Ra) i / 22 ? (s4.9) 
= Re nce ; Rea Ve 
2 a | | a 

2 =? 4(2nf) C Re ge) Cc + (nf) (e Ca)” z Cs Zi , (84.10) 


ajere ijec] on 
Here, e is the equivalent input noise voltage PSD of the CNRS-HEMT C-H and z is their 
equivalent input noise current PSD. e? is the equivalent input noise voltage PSD of the Rear-OPA 
and i is its equivalent input noise current PSD. 

e =2.25 (aV)? /Hz and i? =4(pA)}/Hz in f> 10 kHz [S4R4]. In Eq.(s4.9)-(s4.11), 
(Riv / Ri Jebu is two orders of magnitude smaller than e?. (1+ Ru /R Ù i? / g? is one order of 


magnitude smaller than Z in Eq.(s4.9). (C,, + Cy tC Ra /Ry) @/ gis one order of 


magnitude smaller than C? e4, in Eq.(s4.10). (Cu +Ca + Cy Ru /Ra )( (1+ Ry, /R)i J/g is 
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one order of magnitude smaller than C, e / A>, in Eq.(s4.11). Further ignoring the small quantities 


in Eq.(s4.9)-(s4.11), 


esete/A,, (s4.12) 
R+i+(2af) CEJA, (s4.13) 
ei, = (i.e) = =—j(2nf)C,2/4, . (4.14) 


54.2.2 The equivalent input noise current of the proposed CryoSTM-TIA 
Putting Eq.(s4.12), (s4.13), and (s4.14) into Eq.(s4.5), the equivalent input noise current PSD of 


the proposed CryoSTM-TIA a is obtained as 
2 =i +4k,T/R, +(I/R, +1/R,) IG tefa » + (2nf) (Gare fA) , ($4.15) 


where Cy = CtC; and C= Cat+Ci+Cj. Eq.(s4.5) is Eq.(3.8) in Article. 


[S4R1] A. van der Ziel, Noise in Solid State Devices and Circuits, Wiley-Inter-Science, New York, 
(1986). 

[S4R2] Y.X. Liang, Low-noise large-bandwidth transimpedance amplifier for measuring scanning 
tunneling shot noise spectra in cryogenic STM and its applications, Ultramicroscopy, 234 (2022) 
13466. 

[S4R3] Z.H. Qian, Study on noise models algorithms and matrix descriptions for integrated circuits, J. 
Northeast Norm. Univ. 35 (2003) 41. 

[S4R4] Data sheet of THS4021 OPA, https:/www.ti.com/lit/ds/symlink/ths4021 .pdf.3 
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Supplemental file 5: Estimating the DC tunneling current error for CryoSTM-TIA 


\(,-D/| (ppm) 


R(Q) 


Figure s5-1 For the proposed CryoSTM-TIA, Er = |/;-J|/|J| vs. R (Re[1 MQ, 1 GQ)]) is simulated by TINA-TI 
with V;=1 mV (Black curve) and Vi=5 mV (Red curve). 


For the proposed CryoSTM-TIA, Fig. s5-1 shows the TINA-TI simulation results Er =|Js-/|/|J| vs. R 
(Re[1 MQ, 1 GQ]) with Vi=1 mV (Black curve) and V;=5 mV (Red curve). The simulation results 


show Er<20 ppm, consistent with the calculated results in Article. 
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Supplemental file 6: Estimation of Inv-Amp input offset voltage 


Figure s6-1 DC circuit of the inverting amplifier. 


Table s6-1 Parameters of THS4021 as the Rear-OPA in Post-Amp 


Open loop DC voltage gain aao 97.5 dB 
Input offset voltage drift [vor/°C 15 V/C 
Input bias current Jpp, Jon 3 pA 
Input bias current drift |ipp|/°C, [fon|/°C Estimated as 100 nA/°C 
CMRR 95 dB 
Ac = Aa0/CMRR 2.5 dB 


The Pre-Amp in the left dashed box of Fig.s6-1 is made of CNRS-HEMT C-H. The source of C-H 
is grounded via a variable resistor Rs. The source of C-H is connected to Ri. Ri=R2=1 KQ, Riri =RitRpi, 
Rpi=10 Q, Ri2=R2tRp2, and variable resistor Rp2 can be adjustable from 0 to 20 Q. Ru ~ Riz Ri=1 
KQO. Its input (the gates of C-H) is grounded, and its two outputs O1 and O2 are suspended. 

Adjust Ruı=Rı2. The current that the constant-current source supplies to the amplifier of the Pre- 
Amp is adjusted to /sour=2/as=2 mA. Adjust Rs and Ru, to let C-H at the ideal operating point (las, Vas) = 
(1 mA, 100 mV), and the DC voltages at O1 and O2 are at the same value Ve. The gate-source voltage 
Vgs =-lasRs, and Ve= Vas-Ves. I assume Vgs =-100 mV for the ideal operating point, so Rs= 100 Q, Ru= 


200 Q, and V.=200 mV. 
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The Post-Amp is shown in the right dashed box of Fig.s6-1. Connecting the Pre-Amp to the Post- 
Amp, the inverting amplifier is formed. Adjust Ru and Rp2, to let C-H at the ideal operating point (Jas, 
Vas)=(1mA, 100mV), and the output DC voltage of inverting amplifier is 0. 

For the inverting amplifier, the output voltage drift is vo, the voltage drift at the output O1 is vı, and 
it at O1 is v2. The voltage drift at the output of the constant-current source is vs, and the voltage drift at 
the source of C-H is vs. By the following equations, the input offset voltage drift of the inverting 


amplifier can be obtained. 


(vovi) Rit vov) R= i , (s6.7) 
(vvi )/RL= ga(V1-Vs)-2mVs + ivp , (s6.8) 
(Vo-v2)/Re+ (vi-v2)/RL= v2/Ru +ibn , (s8.9) 
vs/Rs = ga(Vi-Vs) —ZmVs , (s6.10) 
Aa0(V1tVott-V2) + Ac(VitVott tv2)/2 = Vo , (s6.11) 
Vos = Vo/GA0 . (s6.12) 
aao = ~gmRq/(gaRu) . (s6.13) 


where, vos is the input offset voltage drift of the Macro-OPA for 1 °C. The parameters of THS4012 are 
listed in Table s6-1. From the above six equations, vos, vs, vi, and v2 can be solved out as vos 5- 
i(1+gmRs)/(2g¢m) — Vorega/2m, Vs ~iRs/2, and vı ~ v2 ~ i( 1+egmRs)/(2ga). 

it=—8 pA is the current drift of the constant-current source for 1 °C, and |vor|=15 uV is the input 
offset voltage drift of the Rear-OPA for 1 °C. Hence, |vs|~0.4 mV, |vi]~20 mV, and |vos| is about 0.5 
mV. i.e., in Article, the drift of the gate-source voltage of C-H in the inverting amplifier is about 0.4 
mV/°C, the drift of the drain-source voltage of C-H is about 20 mV/°C, and the input offset voltage 
drift of the inverting amplifier is about 0.5 mV/°C. When the temperature fluctuations of the bipolar 
junction transistors (BJTs) in the constant-current source are controlled within 0.0002 °C, the 
fluctuations at the drain of C-H is below 4 uV, and the fluctuation at the source of C-H is below 80 nV, 
so the operating point for C-H is little changed. The power of the constant-current source is about 24 
mW. It is able to control the temperature fluctuations of the constant-current source within 0.0002 °C 


with the well-designed temperature Control System [S6R1]. 


[S6R1] W. Zhang, et al, Design of Precise Temperature Control System for Saturated Absorption 
Frequency Stabilization of DFB Laser, Semiconductor optoelectronics, 41(4) 560-565 (2020). 
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